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We discuss the possibility of tiie non-one-neutron emission ciiannels in the cold fusion reaction '"Zn 
+ -"'Bi to produce the element Z = 113. For this purpose, we calculate the evaporation-residue cross 
sections of one-proton, radiative-capture, and two-neutron emissions relative to the one-neutron emission 
in the reaction '"Zn -I- -"'Bi. To estimate the upper bounds of those quantities, we vary model parameters 
in the calculations, such as the level-density parameter and the height of the fission barrier. We conclude 
that the highest possibility is for the 2n reaction channel, and its upper bounds are 2.4% and at most less 
than 7.9% with unrealistic parameter values, under the actual experimental conditions of [J. Phys. Soc. 
Jpn. 73 (2004) 2593]. 
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1. Introduction 

Heavy-ion fusion reactions provide a powerful tool for syn- 
thesizing superheavy elements (SHE).'-^^ To verify that a new 
element has been created in such reactions, it is necessary to 
identify which evaporation residue is formed when the ex- 
cited fusion product, in competition with nuclear fission, is 
de-excited by particle or gamma-ray emissions to a certain 
stable state. In this paper, we focus on the experiment per- 
formed by Morita et al. to produce the element 113.^'"*' The 
reaction channel of this experiment for the observed two al- 
pha chains was assigned as the one-neutron (In) evaporation 
channel. We investigate here the sensitivity of this assignment 
to different model assumptions. 

The synthesized elements can be identified by measuring 
the alpha-decay chains to known nuclei. However, this is not 
always straightforward in the case of odd-odd evaporation 
residues, such as those in the experiments cited in refs. 3 
and 4. In such alpha-decay chains, one can expect not only 
ground-state-to-ground-state alpha decay but also the involve- 
ment of excited states. Therefore, it is difficult to compare 
directly the observed alpha-decay energy and lifetime with 
the measured values of known products in the decay chain. 
Uncertainties coming from experimental conditions may also 
lead to difficulties in identifications. Morita et al. have re- 
cently performed an experiment on the formation of the nu- 
cleus ^^^Bh,^* which is the grand-grand daughter of the nu- 
cleus ^^^113, and observed its decay modes. From the com- 
parison of the alpha-decay chain of ^^''Bh with that of 113, 
they further confirmed the In reaction channel assignment. 

In addition to these experimental studies, it is desirable to 
calculate the confidence level with which other reaction chan- 
nels can be excluded from the reaction mechanism considera- 
tion. We therefore calculate the probabilities of the competing 
reaction channels relative to the In channel in standard theo- 
retical models. Although the calculation is performed only for 
this special example, this approach can be applied to any other 
experiments for the SHE production based on cold fusion re- 
actions, offering a simple and practical estimator. 

In practice, for the synthesis of the nucleus ^^^113 using the 
209gj(-702jj^ n) reaction,^^ Morita et al. excluded the follow- 



ing possibilities: (i) radiative capture process (zero-neutron 
evaporation channel) leading to the nucleus ^^'^113, (ii) two- 
neutron evaporation channel leading to the nucleus ^''''113, 
and (Hi) one-proton evaporation channel leading to the nu- 
cleus ^^^112. We will quantitatively estimate the probabil- 
ity of these decay channels by calculating the evaporation- 
residue (ER) cross section. 

For such estimations, it is difficult to calculate the non- 
In reaction using the Monte-Carlo simulation for the de- 
excitation process because of its extremely low survival prob- 
ability. In this respect, we employ the "fusion-by-diffusion" 
(FBD) model, ^'^^ which consists of simple analytic expres- 
sions and is well optimized for various cold fusion reactions. 
We extend it to the non-ln reaction. One difficulty in the esti- 
mation is that the de-excitation process is very sensitive to the 
level-density parameter and the height of the fission barrier. 
To reduce these uncertainties, we estimate worst-case scenar- 
ios for relative probabilities by varying substantially the level- 
density parameter and the height of the fission barrier. We 
demonstrate below that the assignment given to the experi- 
ment cited in ref. 3 is highly probable. 

This paper is organized as follows: In §2, we describe de- 
tails of models and parameters for numerical calculations. In 
§3, we show the calculated results and discuss the probability 
of the radiative capture, two-neutron evaporation, and one- 
proton evaporation channels. Finally, we summarize in §4 the 
confidence level with which the evaporation residues of the 
experiments cited in refs. 3 and 4 can be assigned to the nu- 
cleus 278113. 

2. Model and Input Parameters 

In this study, all calculations are based on the (FBD) model 
proposed by Swi^tecki et al.^'''^ and its straightforward ex- 
tensions, because their model and parameters have been well 
optimized for various cold-fusion reactions and reproduced 
well the experimental data.^ The evaporation residue (ER) 
cross section cter is given by 



O'ER = CTe; 



(1) 
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Table I. Fission barrier heights B{ obtained from ref. 23 and one- 
neutron and one-proton separation energies Si„ and Sip estimated with 
FRDMgS.'" 



Nuclide 


Bf (MeV) 


Sin (MeV) 


Sip (MeV) 




6.12 


7.37 


0.48 


278113 


6.06 


6.32 




27^113 


6.31 


7.59 




278112 


5.99 


7.14 





where cTciip is the capture cross section, Pdift is the probabiHty 
that two touching nuclei can reach the compound state, which 
is descried as "diffusion" driven by thermal fluctuations of 
collective degrees of freedoms, and W^uv is the survival prob- 
ability of the excited compound nucleus. 

One characteristic feature of this model is that the 
neutron decay width is calculated with the "transition 
state" method,^' '^"'^^ rather than the traditional "detailed- 
balance" method.'**^ They showed that the neutron evapora- 
tion width evaluated by the detailed-balance method with the 
temperature-dependent shell correction energy is incorrectly 
suppressed near the fission threshold.'''' We apply the transi- 
tion state method to the proton evaporation channel as well. 

For the nuclear mass, the finite-range droplet model 
(FRDM95)'^' and the Hartree-Fock-Bogoliubov model 
(HFB02)^''' have been recommended in the reference in- 
put parameter library (RIPL-2).^'' In the FBD model, the 
Thomas-Fermi model developed by Myers and Swi^tecki 
(MS96)^^' was used. Among those models, we chose 
FRDM95 as the main model in this work, because not only 
the mass excess but also the shell correction energies and the 
fission barrier heights are consistently provided. Although we 
use the parameters optimized for FRDM95 in all calculations, 
we investigate the worst-case scenario at the different exci- 
tation energies estimated with the three mass models, which 
also causes uncertainty in this study, as will be discussed in 
§2.1. 

We calculate the level-density parameters a with eq. (AlO) 
in ref. 7. The uncertainty of the level density with the shell 
correction energy of FRDM95 was extensively investigated in 
a previous study cited in ref. 21. In RIPL-2, it was shown that 
the expected uncertainty for an unknown level-density param- 
eter is about 6.5%. We thus vary the asymptotic level-density 
parameter by up to 10% (see also Fig. 6.4 in ref. 21). 

We use the height of the fission barrier tabulated in ref. 23 
and the one-neutron and one-proton separation energies cal- 
culated with FRDM95. The values used in the calculation are 
tabulated in Table I. We also analyze the neutron separation 
energies calculated with the other two models and observe 
a small difference. For the uncertainty of the fission barrier 
height, the root-mean square deviation of the calculated re- 
sults from the experimental data is shown to be 0.99 MeV in 
refs. 23 and 24; thus, we increase or decrease the height of the 
fission barrier by + 1 MeV. We also simultaneously increase or 
decrease the shell correction energy, because the fission bar- 
rier in the superheavy-mass region is mainly due to the shell 
correction energy at the ground state. 

For calculations of Pdm, we estimate the temperature at the 
saddle point with the fission barrier heights given in Table I. 
All the other calculations are the same as those in ref. 7. Be- 
cause the fission barrier heights differ from the original cal- 
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Fig. 1. (Color online) Evaporation residue cross sections for the In re- 
action channel in various cold fusion reactions calculated with the FBD 
model. The solid and dashed lines respectively indicate the calculated re- 
sults with and without the correction due to the tai'get thickness. The solid 
circle indicates the experimental data obtained from refs. 25-28. 



culations in ref. 7, we need to fine-tune the neck parameter 
s, which is only one adjustable parameter in the FBD model. 
We choose s - 2.1 fm in order to fit the experimental data 
for the cold fusion reactions performed by RIKEN and obtain 
an overall fit for the experimental data, as shown in Fig. 1. 
In the figure, the solid and dashed lines respectively indicate 
the calculated results of the ER cross section with and without 
the target correction due to the target thickness under actual 
experimental conditions (as will be discussed in §2.7). For 
comparison with the non-ln reaction channels, we fix all pa- 
rameters for the In reaction channel. 

2.1 Range of excitation energies 

One large uncertainty in estimations of the excitation en- 
ergy of a formed compound nucleus is due to large differences 
between various mass models. For comparison, we calculate 
the excitation energies of the compound nucleus ^^^113 in the 
209Bj + 702;ji reaction with FRDM95, MS96, and HFB02. For 
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the mass excesses of ^"^^Bi and ^"Zn, we use the experimen- 
tal masses obtained from Audi2003.^^' Since the experimen- 
tal incident beam energy is 261.4 MeV, the resultant excita- 
tion energies are 15.76, 14.11 and 13.61 MeV forFRDM95, 
MS96, and HFB02, respectively. Note that in the experiment, 
the uncertainty of the incident beam energies is +2 MeV due 
to the target thickness.^' We thus vary the excitation energies 
in the range from 1 1.61 to 17.76 MeV in this study. 

2.2 Proton decay width 

We extend the transition state method to the proton decay. 
The particle decay width is then given as 



or 

~ 2^ Jo 



p^\u*-Bp-Vc-e)de, (2) 



where Bp is the binding energy of an evaporated proton, Vc is 
the height of the Coulomb barrier, and pp is the level density 
for the final state. The symbol U* is the effective excitation 
energy, which is defined by U* = E* — A, where E* is the 
excitation energy and A is the pairing energy normaUzed to 
zero for odd-odd nuclei as given explicitly in APPENDIX B 
of ref. 7. The symbol D is l/pi, where pi is the level density 
for the initial state. 

In the calculation, the height of the Coulomb barrier, Vc, 
has a large uncertainty, due to the difficulty in its direct mea- 
surement. In the standard statical model, Vc is estimated from 
the fusion barrier of its inverse reaction. The empirical fusion 
barrier Bfus is often taken as^"' 

1.44(Z- 1) 



fus 



(3) 



1.18(A- 1)1/3 + 3.928 

However, Moretto pointed out that, in a heavy-mass system, 
the height of the Coulomb barrier for a proton emission is 
lower than that of the inverse reaction due to the polarization 
of many protons,^!' and Parker et al. proposed a new empir- 
ical formula extracted from the experimental systematics,^'^-' 
given by Bf^, = 0.106Z - 0.90. For the nucleus ^'^'^113, the 
heights of the Coulomb barriers are 13.87 and 1 1.08 MeV for 
the standard and new empirical formulas, respectively. In this 
paper, we use Vc = 11.08 MeV to discuss the worst-case sce- 
nario. 

2.3 y-Ray decay width 

For the gamma-decay width, we obtain equations and their 
parameters from ref. 33; however, for simpUcity, we ignore 
the angular-momentum dependence and consider only the El 
transition. The gamma-decay width is thus given by 

(4) 



r 

Ty = h I Wy(U*, ey)dey, 

Jo 



where the transition probability Wy for the El ganama-ray 
emission is given by 

^AU\ey) = -^-^^^e^MiSy). (5) 

The function /ei is the gamma- strength function. FoUowing 
ref. 33, /ei was chosen as the following Lorentzian function: 



/ei = 3.31 X lO-^MeV 



-1' 



fyF 



A {El-e^f + e^r^ 



(6) 



where we consider the values of the width F and the resonance 
energy £0 calculated using the droplet model:^"*-* F = 5 MeV 
and £0 is given by 

Eq = 167.23A-'/\l.959 + U.014A-^'^)-^'^MeV. (7) 

2.4 Survival probability for one-proton emission 

We calculate the survival probability for the In reaction 
with the FED model and extend it to the Ip and 2n reactions. 
In ref. 7, the survival probability for the In reaction, 'Wll^\ is 
expressed as the product of the probability for the one-neutron 
evaporation and the survival probabiUty against the second- 
chance fission or the second-chance neutron emission, 
such that we extend it to the survival probability for the Ip 
reaction, W*^, given by 



(8) 



FTot(£*) 

where Fxot is the total decay width defined by Fxot = Fn + Ff -h 
Tp+Ty. By assuming a standard proton evaporation spectrum 
proportional to k exp(-k/T), where k is the neutron kinetic en- 
ergy and T is the temperature of the transition state for proton 
emission, is given by 



^^(2nd) ^|(l + f)exp(-f) K^O, 
^ 1 1 K<0, 



(9) 



where K = U* - Bp - Vc - E^ and E^ is the threshold en- 
ergy for the second-chance fission or second-chance neutron 
emission whichever is lower. 

2.5 Survival probability for two-neutron emission 

For the 2n reaction, we assume that the survival probabil- 
ity is given by the product of F^ "Vrxm^ for the first neutron 
evaporation, F^^^Vr^"* for the second neutron evaporation, 
and the survival probability against the third-chance fission or 
neutron emission. Since the second neutron evaporation width 
depends on the kinetic energy of the first emitted neutron, we 
approximate the probabihty of the second neutron emission 
using the ensemble average in the first neutron emission spec- 
trum. Wsuv^ is thus given as 



r^^°^(£*) /Fr(£:*-Bn-e)„K3rd),^. 



Wluv (£*) 



-(2n) 



Fn^>)Fru--.„-<.» ^ 

rTot(^*)4ot(^*-^n-<e>) 

(11) 



4of(£*)\4o?(^'-'Bn-^) 



-(211). 



- B„ - e)j 
(10) 



where W'"^'"^^ is the probability for the third-chance fission or 
neutron emission, such as that given in eqs. 8 and 9, < • • • > 
is the ensemble average, and e is the kinetic energy of the 
second emitted neutron. The quantity < e > can be obtained 
as < e >= 2T, assuming the Boltzman distribution for the 
neutron evaporation spectrum. 



2.6 Survival probability for radiative-capture reaction 

In an analogy with the calculation of the In reaction, 
we calculate the survival probability Wsu"' for the radiative- 
capture reaction as the product of the probability of gamma 
emission and survival probability against the second-chance 
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Fig. 2. (Color online) Evaporation residue cross sections for one-proton 
and one-neutron emissions versus excitation energy for the compound nu- 
cleus ^''113. The solid and dashed lines respectively indicate the calcu- 
lated results for the In and Ip reaction channels with the standard level- 
density parameter and fission barrier height. The dotted line shows the cal- 
culated result of the Ip reaction channel for the worst-case scenario with 
Qp X 1 . 1 . The dark-gray, black, and gray an'ows indicate the excitation ener- 
gies corresponding to the experimental incident energies calculated using 
FRDM95, MS06, and HFB02, respectively. The light-gray area denotes 
the incident-energy distribution for the worst-case scenario. 



fission or neutron emission after a gamma-ray emission, 



H'< , given by 



We can calculate W<'"* with the probability below the fission 
or neutron emission threshold energy after a gamma-ray emis- 
sion using the gamma-ray spectrum of eq. 5. W^.°"' is thus 
given by 

rU* 

On) J£ih ' ' y y 



rTot(£*) 



(12) 



rU' 

Jo ^y(U*, £y)dey 



(13) 



where E± is the threshold energy for the second-chance fis- 
sion or neutron evaporation, whichever is lower. 

2. 7 Expected value relative to In reaction channel 

To discuss the possibility of the non-ln reaction channel 
under actual experimental conditions, we calculate the ex- 
pected value relative to the In reaction channel, taking into 
account the effect of the target thickness, given by 



P = 



J^^"" a-ER{E*)dE* 
al$iE*)dE* 



where E* and E* are the upper and lower bounds of the 

max mm t'f 

incident-energy distribution due to the target thickness, re- 
spectively. One advantage of this equation is that we can par- 
tially reduce the effects of uncertainties in cTcap and Pdiff de- 
fined in eq. 1 . 

3. Result and Discussion 

3.1 Probability for one-proton evaporation channel 

To discuss the probability for the \p channel relative to 



Fig. 3. (Color online) Evaporation residue cross sections for the radiative- 
capture and one-neutron emission channels versus excitation energies of 
the compound nucleus ^"ll3. The solid and dashed lines respectively in- 
dicate the calculated results for the In and On reaction channels with stan- 
dard parameters. The dotted line shows the results with fly X 1.1 and a 
decrease in fission barrier height. The arrows and the gray area are the 
same as those in Fig. 1 . 



the In channel, we calculate cr'^ and compare it with cr\ 



Figure 2 shows the calculated result versus excitation energy 
for the nucleus ^^^113. The solid and dashed lines indicate 
the results for the In and Ip reaction channels with the stan- 
dard level-density parameter and fission barrier height, re- 
spectively. For the standard parameter set, the excitation en- 
ergy of the worst-case scenario is estimated to be 17.76 MeV 
with FRDM95, and the resultant E* and E* . are 19.76 and 

' max mm 

15.76 MeV, respectively. The incident-energy distribution is 
denoted by the light-gray area in the figure. In this case, we 
determine the expected value of the Ip reaction channel rel- 
ative to the In reaction channel, p(ip/i"), to be 8.51 x 10""^, 
which is a negligibly small probability for one-proton evapo- 
ration. 

To determine the uncertainty originating from the level- 
density parameter and fission barrier height, we calculate cr^^ 
with flp X 1 . 1 and a decrease in fission barrier height of 11 2 
by 1 MeV (We also change consistently the shell correction 
energy at the ground state). The calculated results are indi- 
cated by the dotted line in Fig. 2. Even in this case, the ex- 
pected value of one-proton emission relative to one-neutron 
emission is less than 4.97 x 10"^. The branching ratio of one- 
proton emission is thus negligibly small. 

3.2 Probability for radiative-capture channel 

Here, we discuss the probability of the radiative-capture 
(14) and In reaction channels. We calculate cr^^* and compare it 
with cTg^^ Figure 3 shows the calculated result. The solid and 
dashed lines indicate the calculated results for the In and On 
reaction channels with the standard level-density parameter 
and the fission barrier height, respectively. The excitation en- 
ergy in the worst-case scenario is estimated to be 13.61 MeV 
with HFB02, and the resultant /Smax and E^in are 15.61 and 
11.61 MeV, respectively. For this case, the expected value of 
the On reaction channel relative to the In reaction channel, 

p(0n/ln)js5 10-4. 

We also increase the level-density parameter by 10% to es- 
timate the worst-case scenario. For the radiative-capture re- 
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Fig. 4. (Color online) Evaporation residue cross sections for the 2n and In 
reactions versus excitation energy of the compound nucleus -^'113. The 
solid and dashed lines respectively indicate the calculated results for the 
In and 2n reaction channels with standard parameters. The dashed line 
shows the result obtained with the level-density parameter a„x 1.1 for the 
2n reaction and a decrease in fission barrier height by 1 MeV (We also 
change the corresponding shell correction energy at the ground state). The 
aiTows and the light-gray area are the same as those in Fig. 1. 



action, we do not change the height of the fission barrier of 
^'^113 for the calculation of the second-chance fission, be- 
cause such changes aff'ect the calculated result for the In re- 
action channel. The calculated result is indicated by the dotted 
line in the figure. For this worst-case scenario, f>(0"/iii) jggg 
than 1.70 X 10"^ which is so small that the radiative-capture 
channel can be excluded. 

3.3 Probability for two-neutron evaporation channel 

To study the probability for the two-neutron evaporation 
channel, we calculate cr'^^- To estimate the level-density pa- 
rameter at the saddle point, we use the elongation of the major 
semi axis at the saddle configuration relative to the spherical 
shape, which is 2.5 fm, the recommended value for the de- 
formed nucleus at the ground state in ref. 7. We also deter- 
mine the sensitivity of the calculation to this parameter and 
observe that it is negligibly small. 

We first calculate cr"^^ with the standard parameter set and 
show the results in Fig. 4. The solid and dashed lines respec- 
tively indicate the calculated results for the In and 2n reaction 
channels with standard parameters. In the figure, we see that 
the probability of the 2n reaction cannot be excluded only for 
the excitation energy estimated with FRDM95. The other two 
estimations are well below the threshold energy for the 2n re- 
action, even if we take into account the effect of the target 
thickness. The expected value of the 2n reaction channel rela- 
tive to the In reaction channel, p<2n/in)^ 2.44x 10"^; thus, we 
can safely neglect the possibility of the 2n reaction channel. 

We next discuss the uncertainty originating from the level- 
density parameter and the fission barrier height. We increase 
fln for the 2n reaction channel by 10% and decrease the height 
of the fission barrier of ^'^113 and the shell correction energy 
by 1 MeV. The calculated result is indicated by the dotted 
line in Fig. 4. In the figure, we can see an unrealistic enhance- 
ment of the 2n reaction channel, compared with the fact that in 
all the experiments of the cold fusion reactions, the observed 
peak values of the ER cross section for the In reaction are 



higher than or comparable to that of the 2n reaction. For this 
case, pPn/in) _i 10"^, which is still a fairly small prob- 
abihty for the 2n reaction despite such unrealistic parameter 
values. 

3.4 Probability for total non-ln evaporation channel 

Finally, we calculate the total non-ln branching ratio, 
p(non-in/in)^ using the sum of the probabilities for the proton, 
radiative-capture, and 2n reaction channels. For the standard 
parameter set, we obtain p("on-in/in) - 2.55 x 10"^, indicating 
that the dominant component of the non-ln channels is the 2n 
evaporation channel. 

For the worst-case scenario, we simultaneously use all the 
parameter sets for each channel discussed in the previous 
sections. Although we independently change the model pa- 
rameters for each channel in discussing the worst-case sce- 
nario, their effects on the other channels are very small, be- 
cause we only change the parameters for the different resid- 
ual compound nuclei at the final state and fix those for the 
In evaporation channel in order to reproduce the experimen- 
tal data. With all the parameters for the worst-case scenario, 
we calculate the expected values at the excitation energies 
estimated with FRDM95 and obtain f (ip/i"> = 4.97 x 10-^ 
p(On/in) ^ 9 94 X 10-4, and p(2n/in) ^ 7 92 x 10-^. The to- 
tal non-ln branching ratio is then p("°n-i"/in) = 8.51 x 10"^ 
which is still a fairly small probability. 

4. Conclusions 

We have discussed the possibilities of proton emission, 
radiative-capture, and two-neutron emission channels in the 
de-excitation of the compound nucleus ^^'^113 produced in 
the 2'^'^Bi + ^°Zn reaction reported in ref. 3. To this end, we 
extend the FBD model to those channels. We vary the ratio of 
the level-density parameter to the height of the fission barrier 
to determine its efifect on calculations and determine the upper 
bounds of those channels. 

To discuss the branching ratio, we calculate the expected 
value of the other reaction channel relative to the In reaction 
channel under actual experimental conditions. For the one- 
proton emission channel, we obtain pCp/'") = 8.51 x lO"'* and 
at most 4.97 x 10"^, indicating that the probability of the one- 
proton emission channel can be ignored. For the radiative- 
capture channel, we obtain = 5.14 x lO"'* and at 
most 1.70 X 10"^. We thus consider that the possibility of 
the radiative-capture channel can be excluded. For the two- 
neutron emission channel, we obtain f>(2n/in) _ 2 44 x 10"^ 
and at most 7.91 x 10"^, which is a fairly small probability 
despite the unrealistic parameter values. For the total non-ln 
branching ratio, we obtain pdion-in/i") = 2.55 x 10"^ and at 
most 8.51 X 10"^. We conclude that the 2n reaction is the main 
branch in the non-ln reaction channels. 

For the estimations of the worst-case scenario, the most 
sensitive parameter is the nuclear mass. The development of 
a mass model would largely improve the determination of the 
optimum bombarding energy in synthesizing new elements. 
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